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Abstract 

We study the temporal evolution of a coherent state under the 
action of a parametric oscillator and a nonlinear Kerr-like medium. 

We make use of the interaction picture representation and use an 
exact time evolution operator for the time independent part of the 
Hamiltonian. We approximate the interaction picture Hamiltonian in 
such a way as to make it a member of a Lie algebra. The corresponding 
time evolution operator behaves like a squeezing operator due to the 
temporal dependence of the oscillator’s frequency. We analyze the 
probability amplitude and the auto correlation function for different 
Hamiltonian parameters and we find a very good agreement between 
our approximate results and converged numerical calculations. 
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1 Introduction 

Coherent states were introduced by Schrodinger in 1926 since the early stages 
of quantum mechanics [XJ - These quantum states are characterized by the 
fact that the trajectory of the center of the coherent wave packet evolves in 
time in the same way as a classical harmonic oscillator and its dispersion 
takes the minimum value allowed by Heisenberg’s principle. These facts 
make them especially useful as a link between the classical and the quantum 
theories. 

Much later, in 1963, Glauber introduced the field coherent states, that 
is, coherent states for the electromagnetic field. These states play an im¬ 
portant role in quantum optics[2j 3], The development of lasers made it 
possible to prepare light fields which are very close to the coherent states. 


1 


Glauber showed that the held coherent states can be obtained from any 
one of the three mathematical definitions: (i) as the right hand eigenstates 
of the boson annihilation operator a\a) = a\a) with a a complex num¬ 
ber, (ii) as those states obtained by application of the displacement opera¬ 
tor upon the vacuum state of the harmonic oscillator D(a) |0) = |a) with 
D(a) = exp(aad — a*a ), and (iii) as the quantum states with a minimum 
uncertainty product (Ap)(Ag) = Tl /2 with A q = A p. The coherent states 
obtained from any one of these definitions are identical when one makes use 
of the harmonic oscillator algebra. Subject to a linear interaction a coherent 
state evolves into a new coherent state, that is, they show temporal stability 
[4, 5|. Nieto and Simmons ®nm generalized the notion of coherent states 
for potentials different from the harmonic oscillator with unequally spaced 
energy levels such as the Morse potential and the Poschl-Tellcr potential. 
The resulting states are localized, follow the classical motion and disperse 
as little as possible in time. Gazeau and Klauder [9] proposed a generaliza¬ 
tion for systems with one degree of freedom possessing discrete as well as 
continuous spectra. These states present continuity of labeling, a resolution 
of the identity and temporal stability. Man’ko and collaborators [10] intro¬ 
duced coherent states of an /-deformed algebra as eigenstates of a deformed 
annihilation operator A = af(n ) where f(n) is a function of the number 
operator that specifies the deformation. These states present nonclassical 
properties like squeezing and antibunching. The properties of their even and 
odd combinations have also been studied mm- 

In the presence of a nonlinear interaction, field coherent states evolve 
into non classical states. This can be achieved experimentally by passing a 
coherent state through a Kerr medium resulting in the appearance of dis¬ 
tinguishable macroscopic superpositions of coherent states, the so called cat 

states mm- 

The parametric harmonic oscillator, namely a harmonic oscillator with a 
time dependent frequency, has been studied from several points of view: using 
the method of adiabatic invariants mmmmmm, super symmetric 
quantum mechanics (2T], algebraic methods [22, [23]) an< 4 different approxi¬ 
mation methods m A particularly relevant realization of the parametric 
oscillator is cavity quantum electrodynamics (CQED) where the frequency of 
a given field mode in the cavity can change in time due to the motion of the 
cavity walls or to changes in the dielectric function of the medium [25] . For 
instance, Wineland et. ah [21] analyzed both theoretically and experimen¬ 
tally the loss of coherence caused by fluctuations in the trap parameters and 
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in the amplitude and frequency of the laser beams, heating due to collisions 
with background gas, internal state decoherence due to radiative decay, and 
coupling to spectator levels. 

In this work we consider a nonlinear system corresponding to a single 
mode held propagating in a Kerr-like medium immersed in a cavity with 
a time dependent frequency. In Section [2] we write the Hamiltonian and 
construct its time evolution operator. In section [3] we follow the evolution of 
coherent states under the nonlinear Hamiltonian and analyze some of their 
statistical properties. 


2 Theory 


Consider a parametric harmonic oscillator immersed in a Kerr-like medium. 
Its Hamiltonian is given by: 

H(t) = ^[p 2 + tt 2 (t)q 2 } + H Kerr (1) 


where Q(t) is an explicit time dependent frequency and Hicerr has to do with 
the Kerr-like medium. We can define the usual annihilation, creation and 
number operators as: 


a = 




:(f } 0 q + ip), 


a) = 




- W), 


n = a) a. 


( 2 ) 


where we have set h = 1 and we write the Kerr medium [27] as HKerr — X'h 2 , 
with x a constant proportional to a third-order nonlinear susceptibility x^ 
which is, in general, a small number [28]. To be specific, in what follows we 
will choose Q(t) = Oo[1 + 2 kcos(2Ho^)] [29] with k also a small parameter. 
The Hamiltonian can be written in terms of ad, a and n as: 

Hit) — 9 0 (h + 1/2) + xh 2 + g{t)(a 2 + a t2 + 2h + 1) (3) 


and g(t ) = cos(29 0 f)(l + k cos(29 0 f)). 

The time evolution operator corresponding to the non linear time inde¬ 
pendent part of the Hamiltonian is given by: 

Uq = exp (— iQ 0 t(h + 1/2) — itxn 2 ) (4) 
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and we can write the time dependent Hamiltonian in the interaction picture 
as 

H^t) = g(t) ( e - 2 iQ(h)t a 2 + a t 2 e 2ia(a)f + 2h + l) (5) 

where we have used the fact that f(n)a = af{n — 1 ) and f(n)a1 = a) f{h + 1 ) 
and the effective frequency Q(h) = O 0 + 2 x(l + h) is a function of the number 
operator. 

Notice that the time evolution operator Uq is exact and includes the 
anharmonicity due to the Kerr medium explicitly. The interaction picture 
Hamiltonian given in Eq. [5] is also exact. The operators given in the in¬ 
teraction picture Hamiltonian do close under commutation, however they 
have an explicit time dependence and the Wei-Norman theorem can not be 
applied. Nevertheless the set {a^ 2 ,h + 1/2, a 2 } also forms the basis of a 
Lie algebra (the su(l, 1) algebra) closed under commutation. In order to 
attain a more manageable Hamiltonian that can be written as a linear com¬ 
bination of time independent operators we approximate the exponentials by 
their average value [3UJ, that is, we make the replacement exp[± 2 zfi(h)t] by 
(ao| exp[±2*r2(h)t] |ao) obtaining the approximate interaction picture Hamil¬ 
tonian: 

tfj(t) = g{t) ( e - 2i ^ 0+2x)t d 2 { e - 4ixt{fl) ) + a) 2 e 2i{Qo+2x)t ( e 4ixt{h) ) + 2h + l) ( 6 ) 

where the expectation value is taken with respect to an initial coherent state. 
The resulting approximate Hamiltonian is similar to that of a degenerate 
parametric amplifier m, where a non linear medium is pumped by a strong 
laser inducing the emission and absorption of photon pairs [5T] . 

With this simplification Hj{t) is an element of the Lie algebra with time 
dependent coefficients and the corresponding time evolution operator may 
be written exactly in the product form [32, [33] 

4 4 

71=1 71=1 

with initial conditions a n (t 0 ) = 0 , and we have chosen the ordering Ad = a) 2 . 
Ad = n, A 3 = a 2 and X 4 = 1 . 

The average takes the form: 

e± 2i(n 0 +2x)t (ao | e ±4i xt n| 0;o) = e ±2*(n 0+ 2 x) * exp[ | ao |2 (e ±4ixt _ (g) 
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The complex, time dependent functions a n (t) needed to construct Ujit) are 
obtained from the following set of coupled, nonlinear, ordinary differential 
equations obtained after substitution of Eq. [7] in Schrodinger’s equation 

d\ = —i(fi + 207/2 + 4 aq/ 3 ) ( 9 ) 

«2 = ~i(f2 + 4 « i / 3 ) 

«3 = -*/se 2a2 

d 4 = -i(f 4 + 2ai/ 3 ) 

where the dot means the time derivative. These equations can be solved ei¬ 
ther analytically or numerically, the equation for a 4 (t) being a Riccati equa¬ 
tion and the equations for the other a's(t) can be obtained by integration. 


3 Statistical properties 


3.1 Probability distributions 

Once we have the explicit form for the time evolution operator, we can eval¬ 
uate the temporal evolution of a coherent state |o) by means of 

\a;t) = U 0 (t)U I (t)\a) (10) 


which is given explicitly as: 


a; t) = N a ^ 

l,m 


oe“ 2 inot )\a 4 e 2 * Qot 
m\l\[(l + 2m)!] -1 / 2 


m 

_ e - ix t{l+2m)^ l + 2m) 


( 11 ) 


with 

N a = exp (— itt 0 t/2 T ol 4 T o 2 o 3 — |ck| 2 /2) . 

The probability of finding the k 'th excited state in the distribution at time t 
is given by Pk(or,t) = \(k\a; t)\ 2 . We obtain: 


Pk(ar,t) 


_ [fe/2] 

m =0 


a^{ae a2 ) k ~ 2m 
m\(k — 2m)! 


2 


( 12 ) 


where [7] means the integer part of 7. 
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In Fig. [l] we show the probability distribution as a function of k for three 
different values of time for the case when the Kerr term x = 0 correspond¬ 
ing to a parametric harmonic oscillator. Here the evolution due to the time 
independent part of the Hamiltonian is that of a harmonic oscillator and 
the interaction picture Hamiltonian is a linear combination of the operators 
{d 2 ,cd 2 ,h} so that it’s time evolution operator is similar to a squeezing op¬ 
erator S = exp[^(<^*d 2 — (at 2 )]. 

The Hamiltonian parameters used in this example are k = 0.05, a = 3 + i3 
and times t — 0 (green), t — 2n (blue) and t = 6 tt (red). At the initial time 
the probability distribution is a Poissonian centered at (n) = 18 as corre¬ 
sponds to a usual coherent state. At f = 27 t its width has decreased and it is 
now centered at k — 10 finally at t = its maximum is located at k = 3, its 
width is even smaller and it presents noticeable oscillations after an initial 
bell shape. These oscillations are evidence of the non classicality of the state 
and are due to the quantum interferences in phase space. 



Figure 1: Probability distribution for the evolved coherent state with Hamil¬ 
tonian parameters x — 0, k — 0.05, a = 3 + i3, a 0 = \/l 8 at times t — 0 
(green), t = 2 tt (blue) and t = 6 tt (red). 


3.2 Auto correlation function 

The auto correlation function is defined as the overlap pUJ 

m = («>(0) | <!>(«)} (13) 
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and it takes large values at times whenever the wave packet resembles the 
original one. When the overlap is complete we have a complete revival oth¬ 
erwise we may have fractional revivals when the overlap is a fraction (1 /q) 
of the total probability. The phenomenon of wave packet revivals (complete 
or fractional) has been observed in many experimental situations in atomic 
and molecular systems [ 35 , 35] • 

Using the explicit forms of the time evolution operators Uq and U/, the time 
dependent coherent state can be expanded in terms of the number eigenkets 


n) as given by Eq. 11 and the corresponding auto correlation function is: 


m = e- 


_ „- i sW 2 + „ 4 + «»a J -| 2 p V- (|z| 2 e” ”«<)* (z ,2 ai e ^ iMk+2 ,y 


E 

k,l ,=0 


k\ l\ 


(14) 

When the coefficient of the non linear term y vanishes we deal with a 
parametric oscillator (Fig|2] top) and the auto correlation function for the 
coherent state | z;t) is a periodic decreasing function of time and its explicit 
form can be written as an exponential 


F(t ) x=0 = exp(-ifi 0 t/2 + a 4 + z 2 a 3 - |^| 2 + z * 2 a ^ 2 ^ 01 + \z\ 2 e a2 ~ iaot ). (15) 


As a reference we show in black the temporal evolution for a held coherent 
state | ze lVtot ). 

When the coefficient corresponding to the temporal dependence of the fre¬ 
quency k <1 and that of the nonlinear term is non negligible (and we thus 
have a nonlinear oscillator) we show in Fig. [2] (intermediate) the auto corre¬ 
lation function for a coherent state | z\t). We can see that there are periodic 
fractional and complete revivals with the revival time T rev = Anh/\ E"{no) \ = 
8 n. Notice also that the auto correlation function is symmetric with respect 
to T rev / 2 = 47 T. In Fig. [2] (bottom) we show the case when neither k nor 
y are negligible; we have a parametric nonlinear oscillator. Here the auto 
correlation function shows only fractional revivals. The revival near 8 n is not 
complete and appears at a time slightly earlier than T rev . Notice also that 
the periodicity with respect to t = 47r has been lost. 

In Fig. [3] we plot the absolute value squared of the auto-correlation func¬ 
tion | W(t) | 2 as a function of time for a state \z; t) with Hamiltonian parame¬ 
ters k = 0.25, y = 0.25. In red we present the result obtained using the time 
evolution operator Uj{t) obtained from the approximate interaction picture 
Hamiltoinian Hj{t) and in black the result obtained when the evolution of the 
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a) 



0 2 n 4 n 6 jt Sn 


b) 




Figure 2: Absolute value squared of the auto-correlation function \F(t)\ 2 as 
a function of time for a state \z;t). Hamiltonian parameters, a): k = 0.05, 
X = 0; b): k — 0, x — 0.25; c): k = 0.25, x = 0.25. The black dotted line 
is for a held coherent state |ze* Qot ). In all cases we have z = 2, «o = 2 and 
Hq = 1- 


system is done numerically taking into account the full Hamiltonian given by 
Eq. [5j Notice the almost perfect agreement between the converged numerical 
result and our approximate result. This is an indication of the quality of our 
approximation when dealing with the interaction picture Hamiltonian. We 
stress the fact that the nonlinearity yA 2 has been taken into account exactly 
by means of the time evolution operator Uq. 






























Figure 3: Absolute value squared of the auto-correlation function \F(t)\ 2 as 
a function of time for a state \z]t). Converged numerical results in black and 
our approximate results in red. Hamiltonian parameters: k = 0.25, % = 0.25. 
In all cases we have z = 2, a 0 = 2 and O 0 = 1. 

4 Conclusions 

In this work we have built an approximate time evolution operator for a 
system composed of a parametric oscillator in a nonlinear Kerr-like medium. 
The Hamiltonian is transformed into the interaction picture and as a result 
we obtained a time dependent Hamiltonian that contains the number op¬ 
erator in an exponential. In order to have a more managable Hamiltonian 
we approximate the exponential by its average value taken between a time 
independent coherent state. With this simplification we can write the Hamil¬ 
tonian in the interaction picture as an element of a finite Lie algebra whose 
time evolution operator can be expressed as a product of exponentials. To 
show the quality of our methodology we calculated probability distributions 
and the auto correlation function for a case where neither k nor x are negligi¬ 
ble. We found that the approximate method is consistent with the converged 
numerical results. 
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